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Metallic magnetism is both ancient and modern, occurring in such familiar 
settings as the lodestone in compass needles and the hard drive in comput- 
ers. Surprisingly, a rigorous theoretical basis for metallic ferromagnetism is still 
largely missing^. The Stoner approach perturbatively treats Coulomb interac- 
tions when the latter need to be large^, while the Nagaoka approach incorporates 
thermodynamically negligible holes into a half-filled band^. Here, we show that 
the ferromagnetic order of the Kondo lattice is amenable to an asymptotically 
exact analysis over a range of interaction parameters. In this ferromagnetic 
phase, the conduction electrons and local moments are strongly coupled but the 
Fermi surface does not enclose the latter (i.e., it is "small"). Moreover, non- 
Fermi liquid behavior appears over a range of frequencies and temperatures. 
Our results provide the basis to understand some long-standing puzzles"*'^'^ in 
the ferromagnetic heavy fermion metals, and raises the prospect for a new class 
of ferromagnetic quantum phase transitions. 



A contemporary theme in quantum condensed matter physics concerns competing ground 
states and the accompanying novel excitations^. With a plethora of different phases, mag- 
netic heavy fermion materials should reign supreme as the prototype for competing order. 
So far, most of the theoretical scrutiny has focused on antiferromagnetic heavy fermions^i^. 
Nonetheless, the list of heavy fermion metals which are known to exhibit ferromagnetic or- 
der continues to grow. An early example subjected to extensive studies is CeRu2Ge2 (ref.— 
and references therein). Other ferromagnetic heavy fermion metals include CePtii, CeSix^, 
CeAgSb^i^, and URu2_xRexSi2 at x > 0.3^. More recently discovered materials include 
CeRuPO^ and UIr2Zn2cr^. Finally, systems such as UGe^^ and URhGe^ are particularly 
interesting because they exhibit a superconducting dome as their metallic ferromagnetism 
is tuned toward its border. Some fascinating and general questions have emerged, yet they 
have hardly been addressed theoretically. One central issue concerns the nature of the 
Fermi surface: Is it "large," encompassing both the local moments and conduction electrons 
as in paramagnetic heavy fermion metals^i^, or is it "small," incorporating only conduc- 
tion electrons? Measurements of the de Haas- van Alphen (dHvA) effect have suggested that 
the Fermi surface is small in CeRu2Ge2^'^'^, and have provided evidence for Fermi surface 
reconstruction as a function of pressure in UGe2^^. At the same time, it is traditional to 
consider the heavy fermion ferromagnets as having a large Fermi surface when their rela- 
tionship with unconventional superconductivity is discussed^^i^; an alternative form of the 
Fermi surface in the ordered state could give rise to a new type of superconductivity near 
its phase boundary. All these point to the importance of theoretically understanding the 
ferromagnetic phases of heavy fermion metals, and this will be the focus of the present work. 

We consider the Kondo lattice model in which a periodic array of local moments interact 
with each other and with a conduction-electron band. Kondo lattice systems are normally 
studied in the paramagnetic state, where Kondo screening leads to heavy quasiparticles in 
the single-electron excitation spectrunJ^. The Stoner mean field treatment of these heavy 
quasiparticles may then lead to an itinerant ferromagnet^. With the general limitations 
of the Stoner approach in mind, here we carry out an asymptotically exact analysis of the 
ferromagnetic state. We are able to do so by using a reference point that differs from either 
the Stoner or Nagaoka approach. 

The model contains a lattice of spin-^ local moments (S, for each site i) with a ferro- 
magnetic exchange interaction (/ < 0), a band of conduction electrons (c^^, where K is the 
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wavevector and a the spin index) with a dispersion e^^ and a characteristic bandwidth 
and an on-site antiferromagnetic Kondo exchange interaction {Jk > 0) between the local 
moments and the spin of the conduction electrons. The corresponding Hamiltonian is 



(1) 



The Hamiltonian above is to be contrasted with models for double-exchange ferromagnets 
in the context of, e.g., manganites, where it is the "Kondo" coupling that is ferromagnetic 
due to the Hund's rule. 




}AF exchange 
J >0 



FM exchange 
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FIG. 1: An illustration of the Kondo lattice. Local moments from f-orbitals are in green, and 
are depicted here to be spin down. Spin-up conduction electrons are in red, which have a higher 
probability density than the spin-down conduction electrons in blue. The Hamiltonian for the 
model is given in Eq. ([T|) where a is the spin index and a refers to the three spin directions. Note 
that the Einstein summation convention is used on indices. For simplicity, we assume = 
The characteristic kinetic energy, W, is defined as W = 1/po, where po = Ylx^i^F — e^) is the 
single-particle density of states at the Fermi energy {Ep). Both Ep and the chemical potential, 
/i, scales like W. We use the Shankar notation with K = \K\ measured from the center of the 
Brillouin zone. 

The parameter region we will focus on is Jx ^ |/| ^ W. Here we can use the limit 
Jx = as the reference point, which contains the local moments, representing the f-electrons 
with strong repulsions, and conduction electrons. As illustrated in Fig. [H the local moments 
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order in a ferromagnetic ground state because / < 0, whereas the conduction electrons form 
a Fermi sea with a Fermi surface. A finite but small Jk will couple these two components, 
and its effect is analyzed in terms of a fermion+boson renormalization group (RG) proce- 
dure ^iSi;^. We will use an effective field theory approach, which we outline below and 
describe in detail in the Supplementary Information. Though our analysis will focus on this 
weak Jk regime, the results will be germane to a more extended parameter regime through 
continuity. 

The Heisenberg part of the Hamiltonian, describing the local moments alone, is mapped 
to a continuum field theory^ in the form of a Quantum Nonlinear Sigma Model (QNLcrM). In 
this framework, the local moments are represented by an 0(3) field, m, which is constrained 
non-linearly with a continuum partition function. Combining the local moments with the 
conduction electrons, we reach the total partition function: Z = j T>mV[ip , ip] 5(m^(x, r) — 
1)6^"^, where S = + S'^ + Sk- The action for the conduction electrons, S'^, is standard. 
Defining m"*" = m^^ + irriy and m~ = rrix — iruy, the low energy action for the local moments 
is expressed in terms of a single complex scalar: 

5m ~ - / dud'^q m~^{q,iu){—MQiu + psq'^)m~{—q, —iu) + g f (dm)'^ (2) 
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Here, Mq is the magnetization density, and ps the magnon stiffness constant. The magnon- 
magnon coupling g, schematically written above and more precisely specified in the Sup- 
plementary Information, turns out to be irrelevant in the RG sense when fermions are also 
coupled to the system. Finally, the Kondo coupling can be separated into static and dynamic 
parts. The static order of the local moments induces a splitting of the conduction electron 
band on the order of A ~ J|^(m^) ~ Jf^, which modifies S'^ into the following action for the 
conduction electrons 

S, = J d'^Kde MK, e) (^-le + ^-p + aA^ ^P^iK, e) (3) 

The dynamical part couples the magnons with the conduction electrons, leading to 

The mapping from the microscopic model in Eq. ([1]) to the field theory in ([2])- ([5]) is similar 
to the antiferromagnetic case^S, but differs from the latter in several important ways. One 



simplification is that the translational symmetry is preserved in the ferromagnetic phase. At 
the same time, two comphcations arise. Ferromagnetic order breaks time-reversal symmetry, 
which is manifested in the Zeeman splitting of the spin up and down bands. In addition, the 
effective field theory for a local-moment quantum ferromagnet involves a Berry phase tern>2S 
such that Lorentz invariance is broken, even in the continuum limit; the dynamic exponent, 
connecting uj and q in Eq. ([2]), is z = 2 instead of 1. The effective field theory, comprising 
Eqs. ©-([S]), is subjected to a two-stage RG analysis as detailed in the Supplementary 
Information. 
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FIG. 2: Phase space for the Kondo coupling, a, The spin-splitting of the conduction electron band, 
which kinematically suppresses interband processes associated with the Kondo spin-flip coupling 
to the local-moment magnons. b, The kinematics for the spin-flip Kondo coupling. The low-lying 
excitations of the local-moment system are the magnons which enter the continuum at finite co 
and q. Those of the conduction electrons are expressed in terms of the spin- flip continuum, whose 
Kondo-coupling to the local-moment magnons is cut off below the cutoff energy, cjc « (//bF^)A^, 
and the cutoff momentum, qc ~ Kp-^ — Kpi ~ {Kp/W)A. 

For energies and momenta above their respective cutoffs, Uc ~ {I/W'^)A'^ and qc ~ 
{Kp/W)A, the magnons are coupled to the continuum part of the transverse spin excitations 
of the conduction electrons, see Fig. [21 Here, the Kondo coupling is relevant in the RG sense 
below three dimensions. This implies strong coupling between the conduction electrons and 
the local moments, and both the QNLctM as well as the action for the conduction electrons 
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will be modified. Explicitly, the correction to the quadratic part of the QNLcxM is 



where 7 is a dimensionless constant prefactor. At the same time, the conduction electrons 
acquire the following self-energy: 



where A2 and A3 are dimensionless constants of order unity. Similar forms for the 
self-energies appear in other contexts, notably the gauge-fermion problem and the spin- 
fiuctuation-based quantum critical regime. The formal similarities as well as some of the 
important differences are discussed in the Supplementary Information. 

With these damping corrections incorporated, the effective transverse Kondo coupling, 
J^, becomes marginal in the RG sense in both two and three dimensions. This signals the 
stability of the form of damping for both the magnons and conduction electrons2ii2I. At the 
same time, the effective longitudinal Kondo coupling, J^, as well as the non-linear coupling 
among the magnons, g, are irrelevant in the RG sense. 

The marginal nature of the Kondo coupling in the continuum part of the phase space 
implies that the effective coupling remains small as we scale down to the energy cutoff ~ tUc 
and, correspondingly, the momentum cutoff q ^ Qc- Below these cutoffs, the transverse 
Kondo coupling, which involves spin flips of the conduction electrons, cannot connect two 
points near the up-spin and down-spin Fermi surfaces; see Fig. [2l Although there is no 
gap in the density of states, as far as the spin-flip Kondo coupling is concerned, the system 
behaves as if the lowest energy excitations have been gapped out. The important conclusion, 
then, is that the effective transverse Kondo coupling renormalizes to zero in the zero-energy 
and zero- momentum limit. This establishes the absence of static Kondo screening. Hence, 
the Fermi surface is small, and this is illustrated in Fig. [3^. 

Our result is surprising given that the ratio Jr/^c ~ / (/Jx(^^)^) ^ 1- By contrast, 
the standard Kondo impurity problem with a pseudo-gap of order Ap^ ^ Jr in the con- 
duction electron density of states near the Fermi energy would be Kondo-screened^>^. The 
difference is that, in the latter case, the Kondo coupling renormalizes to stronger values as 
the energy is lowered in the range Ap^ <^ u <^ W; for J^/Apg ^ 1, the renormalized Kondo 
coupling is already large by the time the energy is lowered to ~ Apg. 




(6) 




d = 2 



d = 3 
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The small Fermi surface we have established is to be contrasted with the large Fermi sur- 
face of a ferromagnetic heavy fermion metal in the Stoner treatment, illustrated in Fig. [3)d. 
In the latter case, the local moments become entangled with the conduction electrons as 
a result of the static Kondo screening. Kondo resonances develop and the local moments 
become incorporated into a large Fermi surface. This Fermi surface comes from a Zeeman- 
splitting of an underlying Fermi surface for the paramagnetic phase; the latter is large, as 
seen through a non-perturbative proo^ that relies upon time-reversal invariance. 




FIG. 3: Contrasting the small and large Fermi surfaces. The spin- up electron Fermi surface is 
drawn in red and larger than the spin-down electron Fermi surface in blue. The larger Fermi surface 
has been made slightly transparent to reveal the smaller sheet, a, The local moments are not part 
of the Fermi surface, b, The static Kondo screening has caused the Fermi surface to expand to 
accommodate the Kondo resonances associated with the local moments. 

The region of validity of Eqs. fl6|7|) corresponds to c^c ^ ^ 1-^1 and Qc ^ q ^ 2Kp. 
This range is well-defined, given that A Jj^{m^) < Jx and that we are considering 
Jk <^ \I\ ^ W. In this same energy and, correspondingly, temperature ranges, other 
physical properties also show a non- Fermi liquid behavior. In two dimensions, the specific 
heat coefficient, C/T ~ T^^^'^ and the electrical resistivity p ~ T^/^. In three dimensions, 
C/T ~ log(l/T) and p ~ T^/^. 

Our result of a stable ferromagnetic metal phase with a small Fermi surface provides 



7 



the basis to understand the dHvA-measured^i^>^ Fermi surface of CeRu2Ge2, which is ferro- 
magnetic below Tc = 8 K. Our interpretation rests on a dynamical Kondo screening effect 
that turns increasingly weak at lower energies. This is supported by the observation of the 
collapsing quasielastic peak measured in the inelastic neutron-scattering cross section as the 
temperature is reduced^. It will be very instructive if the Fermi surface of UGe2'^- is further 
clarified and if systematic dHvA measurements are carried out in other ferromagnetic heavy 
fermion metals as well. With future experiments in mind, we note that our conclusion of a 
small Fermi surface also applies to ferrimagnetic order. 

In the parameter regime we have considered, the non-Fermi liquid features are sizable. 
For instance, the non-Fermi liquid contribution to the self-energy [Eq. ([7])] is, at the cutoff 
energy ujc-, larger than the standard Fermi liquid term associated with the interactions among 
the conduction electrons. It remains to be fully established whether the non-Fermi liquid 
terms in the electrical resistivity and specific heat can be readily isolated from contributions 
of other processes. Still, there is at least one family of materials, URu2-xRexSi2 at x > 0.3, 
in which non-Fermi liquid features have been shown to persist deep inside the ferromagnetic 
regime^. Whether this observed feature is indeed a property of the ferromagnetic phase, or 
if it is related to some quantum critical fluctuations or even certain disorder effects, remains 
to be clarified experimentally. We hope that our theory will provide motivation for the 
experimental search of non-Fermi liquid behavior in ferromagnetic heavy fermion metals as 
well. 

Finally, the existence of a ferromagnetic phase with a small Fermi surface raises the 
prospect of a direct quantum phase transition from a Kondo-destroyed ferromagnetic metal 
to a Kondo-screened paramagnetic metal. This, like its antiferromagnetic counterpart^i^i^^, 
in turn raises the possibility of a new type of superconductivity; the underlying quantum 
fluctuations would be associated with not only the development of the ferromagnetic order- 
but also the transformation of a large-to-small Fermi surface. Accessing the quantum phase 
transition requires that our analysis be extended to the regime where the Kondo coupling 
is large compared to the RKKY interaction, and this represents an important direction for 
the future. 

To summarize, we have shown that the ferromagnetic Kondo lattice has a parameter 
range where the Kondo screening is destroyed and the Fermi surface is small. This con- 
clusion is important for heavy fermion physics. It allows us to understand a long-standing 
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puzzle on the Fermi surface, as epitomized by the dHvA measurements in CeRu2Ge2. It also 
sharpenes the analogy with the extensively studied antiferromagnetic heavy fermion metals, 
where the dichotomy between Kondo breakdown and conventional quantum criticality is 
well established. More broadly, the present work has led to one of the very few asymptoti- 
cally exact results for metallic ferromagnetism. Our findings highlight an important lesson, 
namely that correlation effects can lead to qualitatively new properties even for magnetism 
occurring in a metallic environment. This general lesson could very well be relevant to a 
broad array of magnetic systems, including the extensively-debated iron pnictides^^. 
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Supplementary Information for "Metallic Ferromagnetism in the 

Kondo Lattice" 



Seiji J. Yamamoto & Qimiao Si 
Department of Physics and Astronomy, Rice University, Houston, TX 77005, USA 

I. KONDO LATTICE AND FIELD THEORY 

We begin with a microscopic description of heavy fermion metals in terms of the Kondo- 
Heisenberg Hamiltonian. 

^ = E 'A.^K. + ^ E s^s^ + E JKsHJ-fc., (1) 

k (ij) t 

where a labels the three spin components. For simplicity, and without loss of generality, 
we will consider only nearest-neighbor {{ij)) ferromagnetic exchange interaction among the 
local moments, and we will also assume eg = By contrast to the purely itinerant 

magnets, the local moments are independent degrees of freedom to begin with and, on 
their own, would be ferromagnetically ordered (/ < 0). These local moments are also 
antiferromagnetically coupled {Jk > 0) to itinerant conduction electrons. The exchange 
interaction among the local moments includes not only the RKKY interaction generated 
by the conduction electrons in Eq. ([T]), but also the RKKY and superexchange interactions 
from other conduction-electron states as well as the direct exchange. 

Since we are interested in the low energy properties of the ferromagnetic phase of this 
system, we adapt an effective field theory previously used for the pure quantum Heisenberg 
ferromagnet^, but extend it to include fermions. Here, the spin is represented by an 0(3) 
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field, m, which is constrained non-hnearly. 



T)mT)['ip, ^jj] 5(m {x, r) — l)e 

= "^'m + '^Berry + 3'^ + Sk 

^ ~ 2 J ^ ^ ^x^' dxf" 



, dm"'(x, t) 



dr 



"^c = j d'^xdr ip^{x, r) - - fi^ ipa{x, r) 
Sk = Jk j d^xdr sl{x,T)m\x,T) (2) 

where, as usual, = ipia-^'ipia' ■ The topological Berry phase term is crucial to get the 



2 

1-2 



dynamics right . If we define the z-axis as the direction of magnetization, we have Vm x A = 
(0,0, 1) = (m) (note that the curl is in field space, not real space). Thus, in a linearized, 
low-energy theory of spin fluctuations, we have A ^ {—iTiy, rrix, 0). Defining = rrix + imy 
and m~ = rrix — imy we obtain a theory of a single complex scalar 

= <Sfn + ^Berry (3) 

We have now arrived at an effective theory of local moment ferromagnetic magnons coupled 
to fermions with effective coupling constant that for simplicity we also label Jk- The map- 
ping from the microscopic model in equation to the field theory in ([2]) parallels the AF 
caseii^. 



II. SCALING ANALYSIS 



We need to carry out an RG analysis for the field theory above several times, both before 
and after self-energies have been incorporated. To begin, we summarize the pure boson 
problem which has been done previously^. The dimension of the m field is fixed by the 
nonlinear constraint m"" {x , T)m"' {x , t) = 1 which requires [m"(x, r)] = 0. In momentum 
space, this becomes [m°'{q,uj)] = —d — Zf,. Unless indicated otherwise, we will exclusively 
be concerned with field dimensions in momentum space, so the arguments will often be 
dropped: [m] = —d — zi,. As usual for purely bosonic RG, the momenta and energies scale 
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simply as [q] = 1 and [u] = Zb, where Zb = 2 is the dynamical exponent for the boson, which 
is consistent with uj ~ g^. The modulo An ambiguity in the Berry phase dictates [Mq] = d, 
and the scale invariance of Sm establishes [ps] = d + Zb — 2. 

Read and Sachdev were the first to point out that higher order gradient terms may be 
relevant. 

•^m ^ 3 j d'^xdr {d^mad^niadymbdynib - 2dymadymad^mbdymi^ (4) 

Using the scaling scheme described above, this coupling, representing magnon-magnon inter- 
actions, has scaling dimension [g\ = d—2. This indicates that, for d >2, the magnon-magnon 
scattering is relevant. We will see later why this term becomes irrelevant when fermions are 
incorporated. 

In parallel to the pure boson problem, there is a well known procedure for handling 
pure fermion problems within a momentum shell approach^. The essential difference from 
the bosonic RG is that the low energy manifold now consists of an extended surface, the 
Fermi surface, rather than a single point. Scaling should therefore be done with respect to 
this surface, and this may be accomplished by a clever change of coordinates for a simple 
spherical Fermi surface. 

When the action contains both bosons and fermions, the momentum shell RG becomes 
much more complicated. In the special case Zf = 1 and 2;;, = 1, we have extended Shankar's 
approach in a straightforward fashion^. However, such an approach does not work if Zf Zb- 
Another strategy has been proposed by Altshuler, loffe, and Millis^, and we adopt this 
method here. 

Each fermion momentum space integral is decomposed into patches of size A/ in every 
direction so that each patch is locally a fiat space. Scaling is accomplished locally with 
respect to the center of each patch. Momenta are therefore decomposed into components 
parallel and perpendicular {k±) to the vector normal to the Fermi surface at this ref- 
erence point. For example, /g^n^ui^g d'^K = ^patches I-if d'^~^kj_dk\\. Note that some authors 
use an opposite naming convention for components; we follow the notation of Ref.-. A tacit 
assumption of this approach is that the boson does not connect two fermions in different 
patches; this is only justified for forward scattering problems like the one we consider in this 
paper. Bosonic momentum integrals are already constrained to a volume of linear dimension 
Ab, which we assume naturally fits inside the fermionic patch: Af, ~ A/ = A. In this scheme, 
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fermionic and bosonic momenta scale the same way, albeit anisotropically. The assignment 



be different depending on how we incorporate the corrections to the QNLsM and fermion 
actions. The scaling analysis will therefore need to be done anew for each case. 

The introduction of fermions and the choice to use the scaling procedure outlined above 
has an immediate consequence on the way we scale the bosonic action. In the pure boson 
case, we can use [Mq] = d. This comes from the modulo 47r ambiguity of the Berry phase. 
Specifically, since e*^'^'^ = 1, we need iAtiS = i27m, where n is an integer. Therefore S is 
quantized at either an integer or half integer value, and is insensitive to the RG rescaling. 
However, since S = Mq J d'^x = MqL'^, and since [L*^] = —d, we must have [Mq] = d^. But 
the anisotropic scalings we employ in momentum space no longer translate simply to a real 
space analysis. We must therefore abandon these dimension assignments for the pure boson 
problem. Instead, we write the action completely in momentum space and live with the 
understanding that after rescaling, the fields m"'{q,uj) and ip{k,e) no longer represent the 
Fourier transforms of ■m°'{x,T) and tp{x,T). This is nothing new since even in the original 
Wilsonian RG formalism the imposition of a cutoff invalidates the interpretation of 0(g) as 
a true Fourier transform of 0(x). 

A second reason to modify the Read-Sachdev assignments for scaling dimensions in the 
pure boson problem is that the addition of fermions acts as a magnetization sink for the local- 
moment system. Of course, the overall magnetization is still conserved in the ferromagnetic 
phase. Furthermore, we assume there are no valence fluctuations (an implicit assumption 
in writing down the microscopic Kondo-Heisenberg Hamiltonian) so we can still treat the 
local moments as 0(3) spins attached to the lattice, and therefore work with the nonlinear 
field theory. 

The way we fix the scaling dimensions is to define the quadratic action according to: 



of values for [e], [/cy], and [k±] will depend on the form of the quadratic action, and this will 




patches 



(5) 




where, as usual, q± ^ as explained in Ref.-. To ensure that these forms are scale invariant. 
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we make the assignments: 



M = 1 

[k±] = l/zb = 1/2 

= -{3zb + d-l)/{2zb)^-{5 + d)/4 
[m] = -{2zb + d+ l)/{2zb) = -(5 + d)/4 (7) 

This information is used to count dimensions for the Kondo couphng (see figure 1). 

^ Jr j d'^~^q±dq\\dujd'^~^ki_dk\\de ^k+g^i^kdn^q + i^k+qdi^k^rn^ (8) 

The tree-level dimension of the Kondo coupling is now easily found. 



[J|] + 2[d.'^~H:^dkiide] + 2[^] + 



m\ 



- [j±] + - 1 + '^^b _ ^ ^Zb + d-l _ 2zb + d+l 
^ Zb 2zb 2zb 

=^ [4] = {3-d)/{2zb) (9) 

The spin-flip Kondo coupling is relevant in two dimensions, and marginal (at the tree level) 
in three dimensions. Usually, when the Kondo coupling is relevant, we expect the model to 
flow to a strong coupling flxed point where Kondo screening sets in, destroying the magnetic 
order and leading to a paramagnetic phase with a large Fermi surface. This, however, would 
be an incorrect, and inconsistent, conclusion. A proper calculation of the self energies and 
subsequent re-analysis of the scaling dimensions around the appropriate fixed point will 
show that there will never be Kondo screening. 



III. DAMPING CORRECTION TO THE QNLSM AND SCALING 

Our analysis so far has been a little too naive. In particular, it describes the wrong flxed 
point. Note that so far we have not considered the z-component of the Kondo interaction, 

J s^m^, which we refer to as the longitudinal channel. This coupling has two important 
effects. First, it introduces the effect of splitting the spin bands of the conduction electrons. 
Second, when the modified bosonic propagator is inserted into the fermionic self energy we 
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will obtain a non-Fermi liquid form when the Kondo coupling is SU(2) symmetric {J'^ = 
Jr ~ Jk)- What is crucial for this, of course, is that the magnons will remain gapless in 
the presense of the Kondo coupling to the conduction electrons, and we wish to show this 
explicitly. With all this in mind, we present below in some detail the calculation of the 
magnon self-energy, as well as an RG analysis with the modified QNLsM. 





FIG. 1: Interaction vertices 



The first observation is easy to demonstrate. For small fluctuations about the ordered 
state, the longitudinal interaction is approximately / (V^t"^? ~ ~ \m'^m^). where 



we have used the constraint = Vl ^ m+m^ . The "1" comes from the magnetization in 
the z-direction, and leads to a shift in the relative density of the up- and down-spin electrons. 
The reference point for our theory should therefore have a quadratic action for the fermions 
of the form 



d xdr tpcrix, t) { dr /i + crA ipaix, t) 

2m 



(10) 



where A ~ J^{m^) ~ J^. We need to write this in momentum space where it has the effect 



of defining a spin-dependent Fermi wavevector: Kp^- = A/2me(yU -|- crA). Expression ([6]) is 
unchanged except for the new definition of Kpcr- We need to build an effective low-energy 
theory around this fixed point, where there is a gap of size 2 A between the up-spin and 
down-spin bands. This form of the fermionic spectrum is essential to correctly capture the 
damping of magnons via the Kondo interaction. The interaction vertices are represented 
diagrammatically in figure [H while the leading contributions to the self energies are shown 
in figure O The real and imaginary parts of the retarded functions can be calculated exactly. 
For example, the contribution from diagram 11^ is 



qn 



+i^^nsgn(C+,i)e(|c 



q + ir^tSgn(C-,T)0(IC-,Tl - 1) VC',T " 1 
11 



1 



lmUi{q,u) 



qn 



Kp^Qii - |C-,tI)a/i - + ^n0(i - IC+,il)A/i - Ci, ;ii) 
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where we have defined C± a- = ^ '^Tk — ' ^^"^ ^ ^ {+> The region in (cu, g)-space where 
the imaginary part is non-zero is depicted in the main paper. A similar exact expression is 
also available in = 3, but the approximate form is perhaps more useful. The bubble 11^ 
in the regime A ^ ^ vpq ^ /i = Kp/{2me) is approximately: 

n^(g, u;) ^ Jp-J^^ (l + (12) 

where 7^ is a constant prefactor which depends on the spatial dimension, and p^^^ = p'^l 
is the density of states at the Fermi level. In two and three dimensions, the explicit expres- 
sions are Po't^^'' = ^ and Po't^^'' = -^Kpa- The uj/q form of the damping is common to 
a variety of systems; in this case it signifies Landau damping of the magnons with spin 1 
excitations of the fermions. 




Jk J K 

FIG. 2: Self Energies 



To satisfy Goldstone's Theorem, it is necessary for all the pieces of 11 to cancel in such 
a way that the full bosonic propagator emerges in massless form. In the gauge-fermion 
problem, this is a consequence of gauge invariance^. In our case, the cancellation is somewhat 
more subtle. First, note that the diagrams 11'-^ and 11^ are explicitly O(J^). Diagrams 11"^ 
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and n^, however, are both hnear in Jk- This is obvious for 11^, whose calculation is trivial: 

Ren|(g,u;) = - J^rK - n^) 

lmRl{q,uj) = (13) 

The sign difference comes from the fact that there is a four-leg vertex J|- for each spin, but 
the sign of the coupling constant depends on a. The reason why is linear in Jk instead 
of 0{J^) can be seen from a simple calculation at (g = 0, a; = 0), which is non-singular due 

to the different spin indices. After performing the Matsubara sum, 



-^{n,-n,) (14) 

Therefore, when the Kondo coupling is SU(2) symmetric the mass terms cancel and 11"^ -|- 
~ Jjcldl^l/Q and thus x~^{Q:'^^) — q'^ + ldJlrY-> where as usual we have neglected the 
hnear in u term because it is less relevant in the RG sense. This special form of the bosonic 
propagator has emerged in a number of other applications, the most famous example being 
the gauge-fermion problem. We will comment on its consequence a little later. 
With the inclusion of damping, the quadratic action now becomes: 

Sm = J dujd^~^q±dq\\ m"*" ^± + b — ^ m~ (15) 

ded'^~^k±^dk\\ tp^^ (ze — v^k^i — a^k^) i^a (16) 



where and h are simply couplings that control the relative scaling between different 
components of the action. Their dimensions will be chosen to ensure the quadratic action 
is scale invariant. Significantly, in this 2;^ = 3 theory the Berry phase no longer controls 
the dynamics, being instead overwhelmed by the damping term. Physically, this is because 
the magnetization of the local moment system is no longer conserved by itself once it can 
exchange spin flips with the conduction electrons. 
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The scaling analysis now needs to be redone. 

[e] = 1 

= 1/^5 = 1/3 

[a] = 1-1/^^ = 2/3 

[b] = 

[ij] = -i3zk + d-l)/{2zi,) = -{8 + d)/6 

[m] = -{2zb + d + l)/{2zb) = -(7 + d)/Q (17) 

Note that in principle and ipa- could scale differently for different spin projections, but 
because of the way they enter the action, we scale them identically. With these choices, all 
the terms in the quadratic action are scale invariant. The Kondo coupling terms, 

"^K = J d'^'^q^dq\i^dujd'^~^k^dk\\de i)k+q,\'4^Kim~ + ipk+q.ii'k^m^ 
= Jr j d'^^^qi±dqii\dLUid'^^^q2±dq2\\duj2d'^^^k±dk\\de 

X ^k+q,-q2,^'^k,^m^m-^ + 1pk+qi-q2,l'^k,imtl^q2 



18) 



(19) 



are easily analyzed: 



K\ 



[J%] + 2[d'^-^k^dk\\de\ +2[^] + 



m\ 



[4] 



^ d-l + 2zb ^3zb + d-l 2zb + d+l 



[Jk] 



Zb 



2zb 



2zb 



K\ 



{3-d)/i2z) 




(20) 



[Jr] + S[d'^~^k±dkiide] + 2[^] + 2[m] 



[Jr 
{l-d)/z 



^ d-l + 2zb ^3zb + d-l 



Zb 



2zb 



^ 2zb + d + l 
' 2^b 



(21) 



The inclusions oiiojq damping into the quadratic part of the boson action has the effect of 
changing the dynamics from 2;;, = 2 to 2;, = 3, however, there is no change to the dimension 
of the spin-flip Kondo coupling. The longitudinal Kondo coupling is irrelevant for any d> \. 

It turns out that a proper analysis of the fixed point requires insertion of the fermion self 
energy as well^, which we turn to next. 
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IV. ELECTRON SELF ENERGY AND NON-FERMI LIQUID BEHAVIOR 



In addition to the scaling analysis, we have another reason to determine the electron 
self-energy. Aniticipating that the non-Fermi liquid contribution from the Kondo coupling 
to the magnons will be cut off at the energy of order ~ ~ {I /W'^)A'^, we wish to 
ascertain the magnitude of the non-Fermi liquid term at this cutoff scale. This will allow us 
to compare this term with some background Fermi liquid contributions. Since the Kondo 
coupling also occurs in the modified magnon propagator, we present here the calculation of 
the electron self-energy in some detail. 

The leading order contribution to the electron self energy in d = 2 is given by the dressed 
boson, bare fermion and no vertex correction, as depicted in figure [21 

/d'^qduj 
(i? + q, iem + iuJn)x{Q, i^n) 

d'^qdu 1 1 



2 



(27r)3 ie + iuj - ^K+q,a - n(g, i^^) 

^ j2 f (Pqduj 1 I 

^< J (27r)3 ^e + ^uJ- ^x. -^cosO q'- n(g, lu) ^ > 

From the previous section we have the result Il{q,iuJn) ~ ~J'k1^'^- Fo^' the integral over 9 

r^TT 1 27rsgnRe(2) r i 

we use: j^—^ = -^^r^ for any complex z. 

fi> ■ ^ j2 f qdqdu 1 1 , 



(27r)3 Kq/m g2 + M y ie+i^^^^K. _ q 
-^Ji I ^ / d9- ' 



^' / (27r)3 Kq/m qi + M j e+u^ + ico^e 



_ 2 I qdqduj 1 1 2tx sgn(e + u) , . 

Kq 



+ 1 



But in the regime of interest, and with the momentum restricted to ~ Kp-, we have 

[fies to 

qdqduj 1 27r sgn(e -|- u) 



^"xf/nf*^ ^ '^^e self-energy then simplifies to 



(27r)3 Rpq/m q2 + j|^M 



t2 im , p , sgn(e + a;) 
= -Jki — ^ / dq / du^^ ^ (24) 

This integral is a little tricky. First note that the frequency integral should have a cutoff, 
but this is complicated by the presence of the sgn function. It would be incorrect to simply 

21 



shift variables a; — > a; + e. The essential identity we need is: 

f du;f{u)sga{uj + e) ^ 2 [ dujf{u) (25) 

which is only true for even functions: f{Lo) = f{—uj). To see where this comes from, note 
first that for even functions: 

/ dujf{uj) = - dujf{uj) 

J a J —a 

Next, to handle the sgn function we partition the integral into four regions: 

r dujf{uj)sgn{uj + e) = - I 'dwf{uj)+ f du:f{u:)+ f du:f{u:)+ f{uj) 

J-A J-A J-e Jo Je 

where the minus sign is the result of the sgn function. Now we use the identity valid for 
even functions: 

/A /•€ pO pe pA 

dujf{uj)sgn{uj + e) = dwf{uj) - / dwf{uj) + / dujf{uj) + / f{uj) 
-A Ja Je Jo Je 

= 2 /" du;f{u;) 
Jo 

Armed with this identity, the self energy is: 



^2 2im r, , A ,2 7e 
-Jk rn.^9.r^ j2., dqqlog ( 1 + 4 ' 



" {Jher (26) 



Had we used a cutoff on the q-integral, we would have ended up with some unsightly hyper- 
geometric functions whose asymptotic form is the same as above, so it is easier to just set 
the cutoff to infinity straight away. For convenience, wc have so far dropped the stiffness 
{ps) factor in the q^ term of the boson propagator. Reintroducing this factor, and taking 
Ps oc I, we end up with the conduction electron self-energy quoted in the main text, Eq. (6). 
Redoing the calculations for d = 3 is relatively straightforward, although now the integral 

— * 

will be UV divergent. The only difference is that now we set K onto the x-axis since the 4> 
variable is the one that runs from — > 27r. This allows us to use the same identity on the 
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integral that we used in the d = 2 case for the 6 integraL 

q^dq sin OdOdu 1 1 



(27r)4 Kq/mq2 + j2^^\A J ^ l£±i^iii^ _ cos . 
q'^dq sin OdOduj 1 1 /" 1 



^/ (27r)4 Kq/mq2 + j2^^\Aj l±fi^ + , cos . 



-^4 



^2, 



q dq sin 9 dOdto 1 1 27r sgn(e + cj) 

Kq/mq2^j2^^\A ^ , ^ 



Within the regime of interest this simphfies to 



2im 



{2'KfKF-i 

2im 



A'log{l + ^)+Jhelog(l+ 



{2TlfKFl 

So the leading singularity in c? = 3 is: 



47-47loge + Ji7log 



47j 



t + 0{e') (27) 



S oc iJ^eloge (28) 



Again, recovering the stiffness factor leads to the form of the conduction electron self-energy 
presented in the main text, Eq. (6). 

Holstein, Norton, and Pincus were the first to show that the transverse electromagnetic 
field coupling remains unscreened and can in principle lead to non-Fermi liquid behavior-. 
For a real electromagnetic field, the smallness of the fine structure constant suppresses this 
effect to extremely low temperatures. Related non-Fermi liquid form appears in the gauge- 
fermion problem^'^'^. More recently, similar self energies have been found near quantum 
critical points and the nematic fermi fiuidiiii^ii^. The prevalence of this self energy results 
from the generic presence of a massless Zf, = 3 boson coupled to a system with a Fermi 
surface. The problem we have considered here has some important formal differences from 
the gauge-fermion and critical Fermi liquid cases, even in the Zb = 3 continuum regime. One 
difference is in the mechanism by which the boson propagators are gapless. In the gauge- 
fermion problem, gauge invariance guarantees the cancellation of the mass term upon adding 
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the bubble and tadpole diagrams in a large-N calculation of the self energy of the vector 
potential. At the ferromagnetic QCP, the divergence of the correlation length {^^"^ 0) 
leads to gapless quantum critical fluctuations. In our case, it is the SU(2) spin symmetry 
of the Kondo interaction which dictates that the contribution from the longitudinal channel 
exactly cancels that from the transverse channel. A similar effect from the longitudinal 
mode of the ordered itinerant antiferromagnet was recently discussed by^^, and we suspect 
that the cancellation argument we advance here may apply to their well. Another 

feature that is unique to our problem corresponds to the specific non-linear terms [Eq. H] 
that occur here, which come into play in our RG analysis. We have shown that these terms, 
while relevant for the pure Heisenberg problem, become irrelevant when the Kondo coupling 
to the fermions is introduced. 



V. SCALING WITH FULLY DRESSED PROPAGATORS 

Now that we have the expression for the electron self energy we can finally incorporate 
it into the fixed point and redo the scaling analysis. 

^rn = j duj(f~^q^dq\\ m+ + m" (29) 

Sc = j ded'^-^k±dki\ {\ef'^ - vp^ - a^kl) ip^ (30) 

Note that the self energy correction to the fermion in c? = 3 is actually eloge, but for the 
purposes of scaling we can simultaneously treat the cases d = 2 and = 3 by analyzing 
the form e'^^^''. To make every term in the quadratic action scale invariant we make the 
assignments: 

[k±] = l/d 

[h] = 1 

[e] = Zb/d = 3/d 

K] = 1-2/ci 

[t/;] = -{3d + Zb-l)/i2d) = -{3d + 2)/{2d) 

[m] = -{2d + Zb + l)/i2d) = -{2d + A)/{2d) (31) 
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Inserting these dimensions into the Kondo couphng produces: 

[4] = (3 - z,)/{2d) = (32) 
[Jk\ = (3-Z6-d)/d=-l (33) 

In both d — 2 and d = 3, we find that the insertion of the self energies has led to the 
marginality of the transverse Kondo coupling, and the irrelevance of the longitudinal channel. 
This demonstrates that with the correct self energies built into the theory, which references 

the appropriate stable fixed point, there is never any unstable flow of the Kondo coupling. 
The ferromagnetic phase with a small Fermi surface is stable to the Kondo coupling. 
Parenthetically, note that the magnon scattering term scales like: 

9 J {d'-'q±dq\\dujf {q^mf (34) 
3(d - 1 + d + ^b) + 4 - 2(20? + + 1) 



d 

which is always irrelevant. 



l-Zi-2d 

-24^ (35) 



VI. THE EFFECT OF THE CUTOFF 

Below the cutoff, uj < ujc ^ (//H^^)A^ and q < qc ^ {Kf/W)Is,, the transvere Kondo 
coupling becomes irrelevant in the RG sense due to phase space restrictions. The longitudinal 
Kondo coupling, having the scaling dimension (1 — d)/^^, is irrelevant as well. The non- Fermi 
hquid effect will therefore be cut off in this range. 

To ascertain the strength of the non-Fermi liquid contribution, we can compare the con- 
tinuum contribution to the self energy, Eq. (7) of the main text, with the background Fermi 
liquid contribution at the cutoff frequency ujc- Adding a Coulomb interaction u among the 
conduction electrons leads to a Fermi-hquid contribution to the self-energy of the order 
^fl{^) ~ u^Pqc''. In d = 2 we have 

E^^^e ~ a;e) ^ {pMn"^T ~ jf IW^^I' (36) 
^FL{e ~ a;c) = u^pWc ~ {uH^lW')j\, (37) 
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In the parameter range we consider, Jk <^ |/| <^ W , Tinfl{^ ~ <^c) is much larger than 
Sir£,(e ~ ujc). Note that in three dimensions, EjvFL(e ~ (^c) ~ PoJk^c/ I ~ JrI^^i leading 
to a similar conclusion. 
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